Origin of the Sub-diffusive Behavior and Crossover Prom a Sub-diffusive to a 
Super-diffusive Dynamics Near a Biological Surface 
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Diffusion of a tagged particle near a constraining biological surface is examined numerically by 
modeling the surface-water interaction by an effective potential. The effective potential is assumed to 
be given by an asymmetric double well constrained by a repulsive surface towards r = and unbound 
at large distances. The time and space dependent probability distribution P(r, t) of the underlying 
Smoluchowski equation is solved by using Crank-Nicholson method. The mean square displacement 
shows a transition from sub-diffusive (exponent a ~ 0.43) to a super-diffusive (exponent a ~ 1.75) 
behavior with time and ultimately to a diffusive dynamics. The decay of self intermediate scattering 
function {Fa{k, t)) is non-exponential in general and shows a power law behavior at the intermediate 
time. Such features have been observed in several recent computer simulation studies on dynamics 
of water in protein and micellar hydration shell. The present analysis provides a simple microscopic 
explanation for the transition from the sub-diffusivity and super-diffusivity. The super-diffusive 
behavior is due to escape from the well near the surface and the sub-diffusive behavior is due to 
return of quasi-free molecules to form the bound state again, after the initial escape. 

PACS numbers: 



I. INTRODUCTION 

The dynamics of the water molecules in the vicinity of 
a biological surface is found to be different from those 
in the bulk [1 — 10]. The dynamical properties of water 
around a protein surface are found to depend on the dis- 
tance from the biomolecular surface ||l|, ||, ||, ^ . Similar 
observations have been obtained from the computer sim- 
ulations of water near micellar surface ||l0|| . In particular, 
the mean square displacement (MSD) of water molecules 
close to biological surfaces is found to be sublinear with 
time ^, |ri|. These results have been confirmed by 
neutron scattering experiments also ||l^. MSD can in 
general be written as below, 



MSD^ {\Ar{t)f) ^c + a*!" 



(1) 



where |Ar(t)| is the displacement of a molecule at time t. 
For diffusive dynamics, a — 1. For sub-diffusive dynam- 
ics, a < 1 whereas the dynamics is called super-diffusive 
if a > 1. Recent computer simulation studies by Can- 
nistraro et al. showed that the dynamics of water near 
a constrained biological surface is sub-diffusive for short 
time (below 10 ps) with the exponent a varying from 0.75 
to 0.96 depending on the region of water molecules cho- 
sen at different distances from the protein surface ^ . 
Sub-diffusivity has been observed in many other type of 
systems also, such as membrane bound protein Jl^ , poly- 
mer melts jlj, |l^, ^ |l^ , and the transverse motion of 
membrane point due to the bending modes etc. On 
the other hand, enhanced or super-diffusion with an ex- 
ponent 1.5 has been observed in the mean square dis- 



placement of engulfed microsphere in a living eukaryotic 
with a sub-diffusive behavior at short times 
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The occurrence of super-diffusion has been inter- 
preted as the evidence for a generalized Einstein relation, 
whereby the motion of the particles along microtubules 
in a dense network require displacement of the surround- 
ing filaments. The time dependent diffusion coefficient is 
related by time dependent viscosity of the medium|21|. 
The origin of this crossover from sub-diffusive dynam- 
ics at short times (t<10ps) to super-diffusive dynamics 
at intermediate time (10ps< t> 200ps) is not yet fully 
understood. 

Observations for the sub-diffusive dynamics of water 
near the biological surface have often been attributed to 
enhanced stability of the water molecules at the surface 
due to hydrogen bonding to the biological surface |2|, 
p^ , p^ , p5[ . Theoretical models have been developed in 
terms of bound^free equilibration [ p6| . The evidence of 
such exchange was obtained by nuclear overhausser effect 
experiment |M . On the other hand, molecular dynamics 
simulations of water molecules at protein and micellar 
surfaces show the existence of sizable fraction of water 
molecules in the layer with residence time larger than the 
rest [ p8| . The binding energy distribution varies from 
0.5 to 9 kcal/mole. The trajectory of individual water 
molecules clearly shows a dynamic equilibrium, between 
a bound state and a free state [|9[ |o). To the best of 
of our knowledge, no simple explanation has been offered 
for the observed super-diffusive behavior. 

Note that no detailed analytical/numerical solution of 
diffusion in an underlying model potential has ever been 
carried out, even though such a study can throw light 
on the origin of the observed behavior. However, purely 
analytical studies in these type of systems have been lim- 
ited by the fact that the heterogeneous surface faced by 
the water molecules makes a general analytical solution 
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virtually impossible. 

It is shown in recent simulations of a micelle that the 
bound water molecules are energetically stable with a 
potential energy contribution of 9-12 kcal/mole while the 
quasi-free water molecules which surround the bound wa- 
ter are only stable by 5-6 kcal/mole ||3l[. Radial distribu- 
tion function shows that the first shell of water molecules 
shows a peak around 3.5 A while the second shell shows 
a peak around 6-7 A^lj]. The structure of g{r) suggests 
the existence of an effective binding potential which could 
be double well in nature because the effective potential 
is related to radial distribution function by the following 
well-known relation 



(2) 



Recently, Garcia and Hummer reported a study of the 
dynamics of cytochrome c in an aqueous solution by 
molecular dynamics simulation js^] . They observed that 
the mean square displacement of the cytochrome c ex- 
hibits a power law dependence in time with an exponent 
around 0.5 for times shorter than lOOps, and an expo- 
nent of 1.75 for longer times. So they clearly observed a 
cross-over from sub-diffusive in the short time to super- 
diffusive behavior in the longer time. While the short 
time sub-diffusion has been explained in terms of multi- 
basin dynamics, the super-diffusion is explained in terms 
of the enhancement of sampling the configurational space 
due to the concerted motion of atoms. 

In this work, we explain the origin of both the sub- 
and the super-diffusive behavior of the biological water 
by studying the diffusion over simple double well and 
Morse type of potentials. The double well potential rep- 
resents the potential environment of the water close to a 
biomolecule. The first potential energy well corresponds 
to the bound water molecule and the second well corre- 
sponds to the quasi free water molecule. The effective 
potential energy of interaction between the surface and 
the water molecules of course becomes negligible at a few 
layers from the surface and here the water is free and be- 
haves like as in bulk water. A model of free diffusion 
(without any potential) is studied to compare the dy- 
namics of the biological water relative to the bulk water. 
We demonstrate by numerical computation of the prob- 
ability density P(r, t) on the potential energy landscape 
that this model is able to catch the dynamical anomalies 
present in the biological water. The sub-diffusive behav- 
ior is shown to originate not only because of the bind- 
ing to the surface but also due to the bound ^ free dy- 
namic equilibration present among the water molecules 
at the surface. MSD calculated from the distribution 
P(r, t) shows a cross-over from sub-diffusive behavior in 
the short time to a super-diffusive nature in the long time, 
precisely of the type observed in computer simulations 
p2[ . However, we find that a pronounced sub-dijfusivity 
is present only when the exchange between the two wells is 
facile. Self intermediate scattering function Fs{k,t) cal- 
culated from P(r,t) also shows a marked non-exponential 
relaxation. 
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FIG. 1: Potential energy diagram for different models used 
in this study. Model I, II and III are double well potential. 
Model IV is simple Morse potential and Model V is a free 
diffusion model without any potential. 



II. MODEL AND NUMERICAL TECHNIQUE 

Two Morse potentials are combined together to form 
different models of double well potentials. The two wells 
are constructed in such a way so that the separation be- 
tween the wells are tq, 2ro, 3ro in Model I, Model II 
and Model III, respectively. In case of all the double 
well potentials, the depth of the first well is taken as -10 
ksT and the second well is -7 ksT. Simple Morse po- 
tential with an well depth of -10 fc^T is termed as Model 
IV. Simulations have also been done without the pres- 
ence of any potential (free diffusion) which is termed as 
Model V. Figure 1 shows the form variation of the po- 
tentials with distances from a constrained surface. All 
the model potentials reach in the long distance. The 
double well potentials represent the stability of the first 
and second coordination shell of the water molecules sur- 
rounding the biological surface. The models have been 
constructed with different separations of double well to 
monitor the dynamics in the well and also the effect of 
exchange on the dynamics of water. For Model I and 
Model II, the two wells are close and that will facilitate 
the exchange whereas in case of Model III, the exchange 
is hindered. The Morse potential (Model IV), on the 
other hand, has no second well to perform the exchange. 
Dynamics of the probability distribution P(r, t) is gov- 
erned by a simple Smoluchowski equation as given below. 



d_ 

dt 



P{v,t) ^V.D.[V ~ l3¥]P{v,t) 



(3) 



. where, (3 = l/ksT and F is the force obtained from 
the derivative of underlying potential(i<^), F = —\7E. 
Without the angular dependence, the radial part of the 
above equation takes a simple form as given below. 



dP 
'dt 



= D 



d^p 



Qj.2 
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where g = V.F. D is taken to be constant. 

The numerical solution of partial differential equation 
is a non-trivial task. The stability of the solution always 
depends on the equation type, variation of the poten- 
tial and proper discretization of the equation. Here, the 
propagation of the distribution function is solved numer- 
ically by Crank-Nicholson technique . This is an un- 
conditionally stable method for the solution of parabolic 
type of equation. The numerical solution of a parabolic 
differential equation by Crank-Nicholson procedure is de- 
scribed by a simple differential equation as below, 



(5) 



where TL is the total Hamiltonian of the system. The 
formal solution of the above equation is. 



u{t) 



e-^*w(0) 



The above equation can be discretized as below. 



,n+l 



1 - |H 



1 



(6) 



(7) 



where n is the time index and j is the space index. This 
is a half implicit scheme and the above equation can be 
solved by inverting a tridiagonal matrix. Although, the 
procedure is unconditionally stable, proper care should 
be taken with the discretization of time and distance. 
Force has been calculated by using cubic spline method 
where the two Morse potentials are joined to avoid the 
discontinuity. The distance r is scaled by the — 3.5 
A where the first peak of radial distribution function 
appears for the water around a biological surface. The 
time t is scaled by r obtained from the diffusion of water 
{Dy,) in the bulk water at 300 K, r = rl/{ioDw) = 8.5 ps, 
where D„ = 2.5 x cm? / sec. The initial distribution 
of P(r, t) (at t = 0) is taken to be a very steep function 
(nearly delta function) peaked at r = rp. Time step At 
is taken as O.OOOlr and Ar is taken as 0.05ro. Simple 
explicit differentiation scheme also works in this case. 



III. RESULTS AND DISCUSSION 

Figure 2(a) shows the normalized probability distri- 
bution of P(r, t) at time t — 200r for all the model 
potentials. P(j,t) becomes non-Gaussian for the double 
and single well model potentials. Figure 2(6) shows the 
normalized probability distribution for the free diffusion 
model. In case of double well potential r^P{r,t) forms 
bimodal distribution at intermediate time. Mean square 
displacement is calculated from r^P{r, t) using the equa- 
tion below. 



MSD = 
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FIG. 2: (a) Probability distribution r^P{r,t) is plotted for 
short distances. All the models except the free diffusion shows 
the non-Gaussian bimodal distribution, (b) The distribution 
r^P{r,t) is shown for Model V (free diffusion). 



Mean square displacement shows different extent of sub- 
diffusive behavior in the short time depending on the na- 
ture of the model potentials. The cross-over to a super- 
diffusive behavior takes place in the intermediate time 
and finally the dynamics become diffusive. This dynami- 
cal dis-symmetry is depicted in the values of exponent a 
obtained from the fitting of AI SD with Eq. 1 for the dif- 
ferent models in different time windows. Table i shows 
the exponents a obtained for the model potentials in dif- 
ferent time window. 

Table 1 

The signature of the crossover from sub-diffusive to 
super-diffusive dynamics. Exponents from the fitting of 
M SD with Eq. I for different models in different time 
window. 



Model 


Sub-diff 


usive(a ) 


Super- diffusive ( a ) 


0.3 ^ IOt 


0.3 ^ 20t 


20 ^ 200t 


0.3 ^ 600t 


Model I 


0.03 


0.96 


1.76 


1.65 


Model II 


0.58 


0.43 


1.75 


1.64 


Model III 


0.96 


0.97 


1.76 


1.81 


Model IV 


1.03 


1.39 


1.79 


1.68 


Model V 


1.00 


0.96 


1.00 


1.00 



(irr^P(r, t) 



(8) 



It is clear from the values of the exponents in Table 1 that 
in the short time (less than IOt) most of the model po- 
tentials show sub-diffusive behavior. However, the sub- 
diffusivity persists for longer time (upto 20r) in case of 
Model II in which the exchange between the two wells is 
facilitated. All the models show super-diffusive dynam- 
ics in the intermediate time with exponents much higher 
than 1.0. In the very long time (above ~ 200t), diffusive 
dynamics starts setting in. So the exponent decreases in 
the fit for the larger time window. Model V shows sim- 
ple diffusion in the absence of potential with exponent 
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FIG. 3: Sub-diffusive behavior of MSD in the short time for 
all the models. The sub-diffusive behavior is most prominent 
for Model II. Free particle diffusion for Model V is simple 
diffusive. 
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FIG. 4: Super-diffusive time dependence of MSD for all the 
different models. All the models (except Model V) show 
super-diffusive behavior. However, the overall diffusion for 
all the models are much less compared to free particle diffu- 
sion. 



always equal to 1.0. Figure 3 shows the MSD against 
for short time of 20 r. Fitting of Model II with Eq. 1 
with a very low exponent is also shown in the same figure. 
Dynamics of MSD changes to a super-diffusive behavior 
in the longer time (beyond 10 20t) as shown in figure 
4. All the models show a pronounced super-diffusive dy- 
namics although the overall MSD is much less compared 
to the diffusion of the free molecules as given by Model 
V. The exponent a obtained from the fit of the MSD 
are given in Table 1. 

Some nice correlations could be drawn from the change 
in the exponent of diffusion in the sub-diffusive and 
super-diffusive. The more pronounced sub-diffusive dy- 
namics shows a less pronounced super-diffusive behav- 
ior. The slowing down in the dynamics due to the bind- 
ing with the surface is reflected in the self intermediate 



0.96 - 

0.92 - 

0.88 - 

0.84 - 

(i('"o.76 - 



Model I 
Model 11 
Model III 
Model IV 
Model V 



0.6 L 
10 



10 10 

Lot (t/x) 



FIG. 5: Self intermediate scattering function {Fs{k,t)) is plot- 
ted against Logit/r) for all the models at = 2n. Due to 
slow diffusion, all the models (except Model V) show non- 
exponential relaxation. Fs{k, t) for Model II in the short time 
shows the signature of sub-diffusive relaxation (see text). 



scattering function Fs(k, i). Figure 5 shows the non- 
exponential time dependence of the Fs (k, t) for all the 
different models at = 7r/8. For the free diffusion, the 
dynamics of Fs{k,t) is exponential as expected. For all 
the other models, the non-exponentiality originates from 
the very slow relaxation of the system. The pronounced 
sub-diffusive behavior is reflected in the short time dy- 
namics of Fs{k,t) in case of the Model II. Similar sig- 
nature of sub-diffusive dynamics is observed in case the 
second-rank dipole-dipole correlation function of lysoz- 
ime in water from computer simulation jsj]. Mode cou- 
pling theory of glassy liquids provide an explanation for 
the sub-diffusive behavior in terms of (Ar^ (t)) and self in- 
termediate scattering function. Fs{k,t) can be expressed 
from Gaussian approximation as, 



Fs{k, t) = exp 



k^{Ar^{t)) 



exp 



6 



(9) 



(10) 



where A is a numerical constant related to diffusion. 
Clearly for a not equal to unity, Fs{k,t) becomes non- 
exponential. F(^k, t) plotted in figure 4 have been fitted 
to stretched exponential and the exponents obtained are 
considerably less than unity in all the models except for 
the free diffusion which models the bulk water is purely 
exponential which is consistent with the relaxation prop- 
erties observed in bulk water lESl. 



IV. CONCLUSION 

Recent observations regarding the slow orientational 
and translational dynamics of the biological water com- 
pared to the bulk water is mapped to a simple problem of 
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diffusion over a double well potential. The first well cor- 
responds to the first coordination shell where the water 
is strongly hydrogen bonded to the biological surface and 
enthalpically very stable. The second well corresponds to 
the relatively free water which constitute the second coor- 
dination shell. Diffusion over such a potential landscape 
is shown to mimic the recent experimental and simula- 
tion studies showing a cross-over from sub-diffusive to 
super-diffusive dynamics. Sub-diffusive dynamics origi- 
nates not only from the relative stability of the potential 
well, but also from the backward diffusion process from 
the second well to the first well. The slow long time relax- 
ation is reflected from the stretched exponential behavior 
o/Fs(k, t). The present study provides convincing ev- 
idence that the binding to a constrained surface along 
with exchange between the free and bound water could 
be the primary reason for the sub-diffusive dynamics of 
biological water. In view of the results obtained in this 
work, the origin of super-diffusivity can be explained by a 
Heaviside step function representation of MSD as given 



by, 



MSD ^ {\Ar{t)\^) = J deP{e) 



&DT{t - k^j) - a^tH{k 



H{t - fc-/) X (11) 



fb + ^bf 



where P{e) is the probability distribution of the potential 
energy of binding (e). a is the separation between two 
wells (around 3 A), fc^/ and kfb are rate of conversion 
from bound to free and free to bound, respectively. H(t) 
is the Heaviside step function which is equal to unity for 
positive values of its argument and zero otherwise. So if 
the rate of conversion from bound to free is large, then 
the M SD will show a super-diffusive behavior. However, 
the free to bound transition retards diffusion and gives 
rise to sub-diffusive behavior. Note that free to bound 
transition is faster than bound to free because of much 
lower activation barrier in the former case. 
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